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Abstract 

We calculate the quasiparticle effective mass for the electron gas in two and 
three dimensions in the metalhc region. We employ the single particle scatter- 
ing potential coming from the Sjolander-Stott theory and enforce the Friedel 
sum rule by adjusting the effective electron mass in a scattering calculation. In 
3D our effective mass is a monotonically decreasing function of throughout 
the whole metallic domain, as implied by the most recent numerical results. 
In 2D we obtain reasonable agreement with the experimental data, as well 
as with other calculations based on the Fermi liquid theory. We also present 
results of a variety of different treatments for the effective mass in 2D and 3D. 
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1 Introduction 



The evaluation of the Fermi surface parameters has been a cornerstone of the Fermi 
hquid theory since its early years. Precise knowledge of these parameters for an 
electron gas, especially in the metallic domain, is not only a fundamental problem, 
but is also extremely important for physical applications. At present, there is some 
controversy about the three-dimensional results at metallic and over-metallic densi- 
ties. The two-dimensional results are also of significant importance due to the recent 
nonvanishing interest in the 2D physics stimulated by high-Tc superconductivity, the 
fractional quantum Hall effect, as well as the development of 2D electronic devices. 

As is well known from textbooks, quasiparticle excitations can be characterized 
by the renormalization constant Z{k-F) which is related to the residue of the Green's 
function at the Fermi surface, and the quasiparticle effective mass m*. In a simple- 
minded physical picture, 1 — Z{kp) and 1 — m*/m both measure the amount of the 
many-body effects in the electron gas. In this paper we will be concerned with the 
effective (renormalized) electron mass. We will use the effective potential coming 
from the Sjolander-Stott theory [|^, and find the effective electron mass by adjusting 
the effective electron mass in a scattering calculation, so that the Friedel sum rule 
is satisfied. Our approach is "hydrodynamic" in a sense that it does not explicitly 
employ the microscopies of the Fermi liquid, however, it requires the "correct" static 
linear response function as an empirical input . Thus, the Fermi liquid character of 
the electron gas will come in indirectly through the linear response which we use in 
a parameterized form. 

Below we will outline the GW format which is a basis for the majority of cal- 
culations of the Fermi surface parameters. We will compare results of different 
approximation schemes with ours. For both 2D and 3D our results are in a reason- 
able agreement with the most recent calculations based on the Fermi liquid theory, 
as well as experimental data. 
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2 Effective mass in three dimensions 



In this section we will outline the GW format for calculating the self energy, and 
present numerical results for the 3D electron gas. It is well known that in 3D 
the effective mass ratio m*/m is less than unity in the high density limit. The 
high density expansion (r^ ^ 1) was obtained in p|, |^. In the metallic region 
(1 < < 8) there has been some controversy about the behavior of the effective 
mass ratio as a function of the ground state density. The formalism for evaluating 
the self energy part was put forward by Hedin ^ {GW approximation). In a more 
rigorous formulation [Q it can be summarized as follows. The standard starting 
point is the Dyson equation for the Green's function: 

G.(k,.) = (0) ' . . (1) 

with e^^ the unperturbed energy and So-(k, cu) the irreducible self energy. The 
effective mass characterizes the quasiparticle excitation spectrum, 

^ 2m* ' ^ ^ 

and in terms of the self energy is then given as: 

m dT,^(k, to) 



m /m = 1 — 
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k dk 
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(3) 



where m is the unrenormalized (bare) electron mass. The irreducible self energy can 
be approximately expressed in Dyson equation form as 

S.(k, ^) = z I ^ I (q, u;')G^ik - q, ^ - ^') , (4) 

where W function incorporates the many-body effects. In general, it can be approx- 
imated by 

W{q,u;) = -^T{q,u;). (5) 

Here Vg = 47re^/g^ is the (3D) bare Coulomb interaction, e{q, u) is the exact dielectric 
function, and r(g, uj) is the vertex correction H. In the work of Hedin M the random 



phase approximation was adopted by putting F = 1 and using the RPA dielectric 
response in (|): 

W''''\q,uo) = Vk/e''''\q,uo). (6) 

Thus, the W^-function was just an effective RPA screened interaction. The self 
energy was obtained by substituting the VT-function (||) into (H) and using the 
noninteracting Green's function in the right-hand-side of (^. Results of Hedin 
are shown in the Fig. 1 (black triangles). The effective mass ratio assumes its 
minimum at ~ 1, then increases in the metallic region, and becomes greater 
than 1 for > 3. More recent results indicate a totally different behavior. In the 
self-consistent approach of Rietschel and Sham the effective interaction function 
had the same form as in (H), but Eqs. (0) and @ were solved self-consistently. 
Their results (Fig. 1, empty triangles) indicate that the effective mass ratio is a 
monotonically decreasing function of in the whole metallic domain. Similar results 
were obtained by Yasuhara and Ousaka ||^ (Fig. 1, empty boxes) who analyzed the 
Landau interaction function ||^ using analytic fits based on the Monte Carlo results 
of Ceperly and Alder |]TU|. In their work the decrease in the quasiparticle mass 
was due to both spin-parallel and spin-antiparallel parts on the Landau interaction 



function. Another analysis of Landau interaction function based on SSTL was 
carried out in [Q. Finally, a self-consistent scheme of Nakano and Ichimaru with dis- 
function incorporating vertex correction (|^) produced a similar decreasing behavior 
of the effective mass ratio (Fig. 1, circles). Their procedure is rather involved, we 
refer to [|T^ for the details. 

Our vehicle will now be the Sjolander-Stott (SS) theory of the two-component 
plasma As it was recently shown it is essentially a fluid description of the 
electron gas, and does not directly employ the microscopic structure of the Fermi 
liquid theory. It uses the correct linear response of the electron g empirical 
input information. It is well known that this theory is capable of producing reliable 
density profiles around a repulsive impurity. The density profile equation resulting 
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from tlie SS tlieory is 



n 



md,3D 



In 



k + q 
k — q 
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(7) 



Tliis gives tlie induced electron density n™'^ around an impurity of charge Z. Here 
f^^{q) is the induced electron density in the linear response approximation: 



3D/ 



k^q^ 



The static linear response x^^l?) is an input information and in our calculations was 
taken from the parameterization in []TB|. fc, q are in units of the Fermi radius kp = 
(STT^n)^, n being the homogeneous ground state density. This profile relation has 
the correct high-g dependence which is responsible for satisfying the cusp condition 
at the Coulomb source: 

= (1 + ^ ^ 



n. 
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The major shortcoming of the profile relation (j^) in the domain of repulsive impu- 
rities is overscreening by a hole. If the (repulsive) impurity charge is big enough the 
total electron density at the location of the impurity goes negative. For most of the 
possible physical applications this feature has an insignificant effect ||l6l because the 
region of the non-physical behavior is very small and one can simply put the total 
electron density to in this region. However, to be on the safe side, we restrict 
ourselves to the range of small enough impurity charges 

- 0.3 < < -0.03 (9) 

where the overscreening by a hole does not occur and one obtains a reliable density 
profile. We extract the effective single particle scattering potential using the SS 
theory 16|. With this scattering potential we check the Friedel sum rule. 

Friedel sum rule is a condition on the difference of the trace of the logarithm of 
the single particle scattering matrix between the top and the bottom of the band 



17| , p!8| . In our case of the jellium model of electron gas it takes a familiar form 



2 °° 

z = - y^'^M^p) ■ 

^ 1=0 



(10) 



where the factor rii accounts for the angular degeneracy: 

f 2/ + 1 for 3D 
ni = < 

[ 2 - 6ifi for 2D 

The factor of 2 comes from the spin degeneracy, 6i{kp) are the partial wave scatter- 
ing phase shifts at the Fermi momentum. This sum rule has been routinely used to 
adjust free parameters of the effective potential in a self-consistent fashion. In our 



case we take the effective potential from the SS theory [0, |T6|. With this effective 
potential we run a scattering calculation at the Fermi surface, and, having obtained 
the phase shifts Sii^kp), check the sum rule ([T0|). The particles that are scattered at 
the Fermi surface are quasiparticles, not bare electrons. Their mass comes explicitly 
into the scattering calculation. We adjust this effective electron mass in the scatter- 
ing calculation until (0) is satisfied within 0.01% accuracy. This provides us with 
the value for the effective mass. The procedure above is repeated for several values 
of Z from the region (|^) in order to insure that the results are independent of the 
impurity charge within a reasonable range 

Our data for 3D are plotted on Fig. 1 (black circles). They clearly indicate that 
the effective mass ratio is a monotonously decreasing function of the Seitz radius 
throughout the whole metallic domain. In the high density limit they (as well as 
all the other data) converge to the limiting behavior 0, ^. The agreement of our 
results with the most recent 3D calculations suggests that the same procedure can 
be tried in 2D since, in principle, the hydrodynamic model does not distinguish 
between dimensions as long as the correct linear response is employed. Now, we will 
consider the situation in two dimensions. 



3 Effective mass in two dimensions 

Two dimensional calculations based on the many-body formalism have been carried 
out within the same format of GW approximation. In the work of Jang and Min 
T^] the 14^-function (|^) is defined so that the vertex correction F = 1. The dielectric 



function is expressed in a standard way in terms of the local field correction Q{q, u): 

e{q,uj) = l-— — ^ — -. 11 

The notation here is as before, but refers to the 2D quantities: Vq = 27r'^/q, and 



the noninteracting response xoiQ,^) is as in |20|. A further approximation is con- 
ventionally made which consists of replacing the dynamical local field correction by 
a frequency independent one: Q{q,uj) = Q{q). Jang and Min employed different 
parameterizations for the Q (q) together with the noninteracting Green's function in 
(^). Their approach was not self-consistent. The employed approximations were: 
the RPA with Q = 0; the Hubbard approximation (HA) adopted to the 2D by 



Jonson with 



where kp = y2jm is the 2D Fermi radius; and the modified Hubbard approximation 
(MHA) [ID with 



^™(?) = ^\/^^ + fc| + ^^F, (13) 

where the Thomas- Fermi momentum is given by A;tf = 27™eVeF = v^r.A^p. The 2D 
Seitz radius is given by 7r(rsao)^ = \jn. The effective mass ratios produced by these 
approximations are plotted in Fig. 2: Hubbard approximation (empty triangles), 
modified Hubbard approximation (empty boxes), and RPA (empty circles). 

We proceed with the scheme we developed for the 3D case, but adapted to 2D. 
The profile relation coming from the Sjolander-Stott theory takes the form 

where the notation is the same as in (0), but refers to 2D. $(g. A;) results from the 
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angular integration 23 



$(g,fc) = ^ 

TT 
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with ttg^k = TT— — 5 



Here, K{x) and 2-^1(2;) are the complete elliptic integral of the first kind and the 
hypergeometric function, respectively. The induced density in the linear response 
approximation P^{q) is taken from the parameterization |16| of the numerical results 
of Neilson et al. |2^. Just like (|^), {lA) has the correct high-g dependence which is 



responsible for satisfying the cusp condition at the Coulomb source: 
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In order to avoid dealing with the overscreening by a hole we take the repulsive 
impurity charge to be small enough 



-0.1 < < -0.01 (17) 

As in the previous section, we extract the effective potential using the SS theory 
and run a 2D scattering calculation at the Fermi momentum adjusting the effective 



electron mass until the Friedel sum rule (10) is satisfied. Just as in 3D, we repeat the 



calculation for different values of Z within the range ([17|). Our results for the effective 
mass ratio are shown on Fig. 2. They fit roughly in between the results of the HA 
and the MHA. The HA is known to take into account the exchange interaction and 
neglect correlations, while the MHA is an attempt to incorporate correlations. To 
compare with the experimental data we plot the magnetoconductivity measurements 
of Smith and Stiles |2^ (Fig- 2, crosses). The agreement of our results with the 



experimental values is rather reasonable. We also plot the high density expansion 
(Fig. 2, dashed line) obtained by Isihara and Toyoda |26| from the specific heat 
calculation: 

mVm = 1 + 0.043r, . (18) 

This result was obtained by including first- and second-order exchange, as well as 
the ring diagram contributions for small but finite temperatures. 

And last, a few words about our numerical procedure would be appropriate. In 
the solution of and (0) we could use (|) and (p!6D, respectively, to approximate 
the high-g behavior, and after solving in the remaining domain match the solutions. 



This procedure is, however, very tedious and ineffective. We found that the solution 
is not affected if we solve the integral equation on the whole semiinfinite domain 
[0, oo) using the Gauss-rational rule for discretizing the integral. Of course, 2D 
calculations require care because of weaker convergence of integrals. In the scattering 
calculation we used 12 partial waves. 

4 Conclusion 

We have presented our results for the electron effective mass ratio using the effective 
potential coming from the SS theory for two and three dimensions. The calculations 
were carried out for small repulsive impurities where the SS works. The effective 
mass was extracted by enforcing the Friedel sum rule in the scattering calculation 
for the effective potential. We compared our results with the most recent numerical, 
as well as experimental data. 

As for the 3D results (Fig. 1), there has been some uncertainty in the behavior 
of m* /m in the metallic domain. Earlier results predicted that this ratio should 
increase with r^, while more recent results indicate that it is a decreasing function 
of Ts- In our treatment the effective mass ratio is a monotonically decreasing function 
of the Seitz radius in the whole metallic domain. 

In 2D (Fig. 2) we present our results together with the results based on the 
GH^-approximation. Our treatment seems to fall "in between'' the Hubbard and 
the modified Hubbard approximations. Our results give reasonable agreement with 
the experimental data. 

In conclusion, it is interesting to relate the 2D and 3D results between themselves. 
It is well known that the exchange interaction diminishes the effective mass, while 
correlations shift it in the opposite direction. The quasiparticlc effective mass is a 
result of an interplay between the exchange and correlation contributions. It comes 
as no surprise that because correlations are stronger in 2D than in 3D, generally, 
the effective mass is greater in 2D than in 3D for the same value of r^. 
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Figure Captions 



Fig. 1. The effective mass ratio m* /m vs. the Seitz radius for the 3D electron gas. 
Resuhs of Hedin (black triangles), Nakano and Ichimaru (empty circles), 
Yasuhara and Ousaka (empty boxes), present (dark circles), and Rietschel and 
Sham (empty triangles). The dashed line is the high density limit f^, ||. Our 
results indicate that the effective mass is a monotonically decreasing function of 
throughout the whole metallic domain. 



Fig. 2. The same as in Fig. 1. Results of Jang and Min [19| in the Hubbard 
approximation (empty triangles), modified Hubbard approximation (empty boxes), 
RPA (empty circles); present theory (black circles); experimental data of Smith and 
Stiles (crosses). The dashed line is the high density result (|l^) of Isihara and 
Toyoda [BH . 
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